Abstract. BeppoSAX observed the Seyfert 1 galaxy 
NGC 3783 for - 5 days from 1998 June 6 to 1998 June 10. 
The average flux during the observation was F2_iofc e y ~ 
6 x 1CP 11 erg cm~ 2 s _1 . The long exposure provided a 
very high quality spectrum that exhibits all the typical 
features of Seyfert Is. The high energy cut-off of the in- 
trinsic continuum, at E c = 340±fQ° keV, was detected for 
the first time in this object. During the observation the 
source showed modest flux variations (~ 20% in 2-10 keV) 
' and spectral variability that can be completely explained 
f^i . by a change of the intrinsic slope Ar ^0.1. NGC 3783 
' is one of the few Seyfert Is that exhibits a soft emission 
, component in the BeppoSAX spectrum. This soft excess is 
5h ' well reproduced by a black body component with temper- 
^ . ature kT~0.2 keV. A warm absorber is also required to fit 
^ ' the soft X-ray spectrum. Both reprocessing features, the 
, Compton reflection hump and the Fe emission line, were 
£S) ' detected in our observation. In addition to a narrow iron 
line (as observed by Chandra) a broad component with 
a = 0.72±J;|f keV and EW=115±76 eV, is required to fit 
, the line. The observed Compton hump and the Fe line are 
^-f. ' well reproduced with the ionized (lg£ = 2.73±q;q|) disc 
ff~) ■ reflection model of Ross & Fabian. In this case the soft 
C*") ' excess observed below 2 keV can be accounted for by the 
, disc emissivity and no additional soft X-rays component 
' is required to fit the spectrum. 
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^p^, Seyfert; X-rays: galaxies 

o 



A&A manuscript no. 

(will be inserted by hand later) 

Your thesaurus codes are: 

missing; you have not inserted them 



The BeppoSAX broad-band spectrum and variability of the 
Seyfert 1 NGC 3783 

A. De Rosa 1 , L. Piro 1 , F. Fiore 2 , P. Grandi 1 , L. Maraschi 4 , G. Matt 3 , F. Nicastro 5 , and P.O. Petrucci 4 

1 Istituto di Astrofisica Spaziale, C.N.R.,Via Fosso del Cavaliere, Roma, Italy 

2 Osservatorio Astronomico di Roma, Monteporzio Catone, Italy 

3 Dipartimento di Fisica, Universita degli Studi "Roma Tre", Via della Vasca Navale 84, 1-00146 Roma, Italy 

4 Osservatorio Astronomico di Brera, Milano, Italy 

5 Harvard-Smithsonian Center of Astrophysics, Cambridge MA 02138 USA 
Received ; Accepted 



ASTRONOMY 

AND 
ASTROPHYSICS 



1. Introduction 

Intrinsic spectral variability has been observed in sev- 



eral Seyfert 1 galaxies: NGC 5548 ( [Nicastro et al. 200C ), 
NGC 4151 (Piro et al. 2002), IC 4329A ( pone et al. 2000| ), 
NGC 7469 (Nandra et al. 2000), MCG-6-30-15 (Vaughan 
& Edelson 2001). The usual behaviour is a steepen- 
ing of the primary power law when the source bright- 
ens. The favoured model for the X-ray emission ob- 
served in AGN comprises a hot, optically-thin plasma 
(the corona) which up-scatters to X-ray energies the soft 
photons produced in a cold, optically-thick accretion disc 



(Haardt & Maraschi 1991). In addition the optically thick 
disc reprocesses and re-emits part of the Comptonized 
flux producing the characteristic reprocessing features 
(the Compton reflection hump and the Kq Fe line). In 
a disc+corona emission model the intrinsic variability has 
to be associated with a change of some properties either 
of the emitting coro na (temperature kT e , optical depth 



r or geometry, see Haardt, Maraschi fe Ghisellini 1997 ) 



or the accretion disc (flares) with the X-ray flux. A de- 
tailed Comptonization model applied to the BeppoSAX 
long look of NGC 5548 (Petrucci et al. 2000) showed that 
the steepening of the spectrum during a flare in the cen- 
tral part of the observation was probably associated with 
an increase of the soft photon luminosity, the hard photon 
luminosity remaining constant. 

The presence of a broad iron line in the X-ray 
spectra of Seyfert galaxies (Fabian et al. 2000) is sug- 



gested by ASCA ( 


Mandra et al. 1997, 


Tanaka et al. 1995 


Yaqoob et al. 2005 


), BeppoSAX (Guainazzi et al. 1999) 



and XMM (Wilms et al. 2001) observations. Neverthe- 
less a narrow iron line component is clearly detected 



in the Chandra (Kaspi et al. 2001, Yaqoob et al. 2001) 
and XMM observations (Reeves et al. 2001, Pounds 
et al. 2001). This narrow component could be attributed 
either to an optically thin medium (like the BLR) , and in 
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this case no Compton reflection should be associated to 
it, or to an optically thick medium with N# > 10 24 cm -2 
(a s the molec ular torus like that observed in the Sy2 
s, Matt 2001a ), whi ch should contribute to the observed 
Compton reflection ( Ghisellini et al. 1994[ ). At least in the 



Alessandra De Rosa: 



case of NGC 4051 (Guainazzi et al. 1998) it is possible 
to associate the narrow iron line component to a distant 
Compton-thick matter. 

The Broad band spectral variability in bright Seyfert 1 
galaxies program, carried out by BeppoSAX in the last 
3 years, exploits the unique capability of BeppoSAX to 
provide spectral informations up to 200 keV. The primary 
goals of this program are to investigate the correlation of 
the intrinsic parameters (r and E c ) with the luminosity 
and to characterize the diverse components of Seyfert 1 
spectra from 0.1 to 200 keV, by studying their response 
to changes of the intrinsic luminosity. 

The bright Seyfert 1 galaxy NGC 3783 was selected 
as a target in the cited BeppoSAX program because: a) 
it is bright enough that the spectrum up to about 100 
keV (F(2 - IGkeV) > 2 x 10~ n erg cm" 2 s" 1 could 
be measured; b) it is characterized by variations of a 
factor of two or more on the time scale of days. This time 
scale is appropriate to check possible correlations in the 
variations of the spectral components. The extensively 
studied X-ray spectrum of this object is characterized by 
narrow and broad features in addition to the intrinsic 
continuum. Deep Oxygen absorption edges (between 0.5 
and 1 keV) we re detected by ROSAT ([Turner et al. 1993 1) 
and ASCA ( peorge et al. 19984 |George et al. 1998b| 
which also observed a variable warm absorbing gas along 
the line of sight. An OVII emission line, likely originating 
in the same warm gas, was also detected in the first 
ASCA observation (in 1993). The Chandra high energy 
resolution observation ( [Kaspi et al. 200C ) identified, in 
addition to many narrow absorption lines, several weak 
emission lines mainly from O and Ne, suggesting the 
presence of two absorption components (outflowing by 
the central source) , that differ by an order of magnitude 
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in their ionization parameters. A Compton reflection 
hump and a Fe Ka emission line were detected by 
GIN GA ([Nandra fc Pounds 1994Q . A Chandra observa- 
tion (Kaspi ct al. 2001) detected a narrow and unresolved 



the Galactic hydrogen column density (Nfj 



gal _ 



9.6 x 

Murphy et al. 1996[). The fit is very bad 



iron line, with an upper limit to the FWHM of 3250 km 
s _1 . These data could not constrain any broad iron line, 
though they are consistent with the broad line found in 
previous ASCA observations. 

In this paper we present the analysis of a BeppoSAX 
observation of NGC 3783. In section ^ we describe the ob- 
servation and the data reduction. In section ^| we report 
the spectral analysis of the data integrated over the to- 
tal observation while in section ^ the spectral variability 
analysis is presented. The results are discussed in section 
|^. A summary is given in section [(]. 



2. Observations and data analysis 



NGC 3783 



observed 



by 



Piro, Scarsi fe Butler 1995) 



BeppoSAX 
for ~ 



( [Boella et al. 1997 
5 days (t exp ~ 154 ks) from 1998 June 6 to 1998 June 10. 
The LECS, MECS and PDS data reduction followed th e 
standard procedure ( Fiore, Guainazzi fc Grandi 1999[ ). 
The PDS spectra were filtered with fixed rise time. 
We extracted the spectrum within a circular region 
centered on the source with a radius of 4' and 6' for 
MECS and LECS respectively. The background was 
extracted from event files of source-free regions ("blank 
fields"). Spectral models were fitted to the data using 
the XSPEC package. All quoted uncertainties correspond 
to 90% confidence interval for one interesting parameter 
(Ax 2 =2.71). The models included a normalization factor 
for each instrument to take into account possible mis- 
calibrations. We allowed the PDS/MECS normalization 
to vary between 0.77 and 0.95, while the LECS/MECS 
normalization ratio was running between 0.7 and 1 
( Fiore, Guainazzi fc Grandi 1999] ). The journal of the 



BeppoSAX observation of NGC 3783 is shown in Table g. 

In Figureg we show the LECS (0.2-3 keV), MECS (2- 
10 keV) and PDS (13-200 keV) light curves. The source 
showed rapid variability both in the soft (0.2-3 keV) and 
medium (2-10 keV) energy ranges, with the larger ampli- 
tude at low energies. The 2-10 keV and 0.2-3 keV flux 
variations were up to a factor 1.6 and 2, respectively, on 
a time scale of about half a day. 

3. The broad band spectral analysis: 0.1-200 keV 

In this section we present the spectral analysis of the 
data integrated over the total observation of NGC 3783 
(t exp ~ 154 ks). We fitted the LECS, MECS and PDS 
data in the energy range 0.15-3 keV, 1.5-10 keV and 13- 
200 keV respectively. In Figure we show the LECS, 
MECS and PDS spectra and the data/model ratio af- 
ter fitting with (model A) a power law absorbed by 



10^ u cm 

(x 2 /dof=818.6/114, with r=1.58±0.06). Several devia- 
tions of the data from the model are apparent: an excess 
of counts at ~ 0.6 keV, a deficit below 2 keV, an emission 
line at ~ 6.4 keV and another deficit above 100 keV. 

We then fitted the data with a model (B) includ- 
ing the known spectral components already detected 
in NGC 3783: a warm absorber ( Turner et al. 1993 ), 
a Compton reflection hump and an iron emis- 
s ion line (modeled wit h a gaussian) at ~ 6.4 keV 
(Nandra & Pounds 1994). The warm gas, assumed to be 
in photoionization equilibrium, is characterized by two 
parameters: the ionization parameter U, defined as the 
ratio of the number density of ionizing photons to hy- 
drogen nuclei, and the hydrogen column density Nff. 
The description and discussion of the photoionization 
code we employed will be given in Sect. [T^. To fit the 
Compton hump we use here a cold disc reflection model 



(PEXRAV in XSPEC, Magdziarz fc Zdziarski 1995). The 



reprocessed spectrum is characterized by the reflection 
fraction R = AVL/2n i.e. the solid angle subtended by 
the reflecting medium. The inclination angle between the 
line of sight and the reflector is assumed i=30°. For the 
moment neither ionization nor smearing effects (gravita- 
tionals and kinematic) are taken into account to fit the re- 
processed features. We will discuss these cases in Sect. |5 . 3| . 
All the line parameters (energy, intrinsic width and inten- 
sity) are free to vary. The effects of including a narrow line 
component - as observed by Chandra ( Kaspi et al. 2001 ) 
- will be discussed in Sect. 5.3. The result of the fit with 
the model B is still unacceptable (x 2 /dof=189. 6/108. See 
Table |J). The ratio data/model B is shown in the upper 
plot of the left panel of Figure ||. This plot shows still 
some deviations of the model from the data: a deficit of 
counts above 100 keV and an excess at energy less than 1 
keV. We therefore added to the primary power law in the 
model B an exponential cut-off (Model C) . The statistical 
significance of this new spectral component is larger than 
99% (following a F-test probability for the addition of 1 
interesting parameter we find F=8.14, see Table |J). 

The excess at low energies may correspond to OVII 
emission line. As su ggested by the ASCA observation 
( George et al. 1998a ), we fitted this excess with a gaus- 
sian emission line. The model is indicated with D and the 
best fit parameters with the statistical significance of the 
new spectral component is shown in the second column 
in Table |[ The ratio data/model D is shown in the lower 
plot of the left panel in Figure ||. 

Despite the good x 2 , there is a problem with model D 
regarding the origin of the emission line at low energy. We 
measure EW oa:!/ g e „=240±27 eV (see Table |). Can this 
spectral feature be produced in the same gas which pro- 
duces the O edges at E < 1 keV (the warm absorber)? 
We suppose the warm absorber to be ionized by the cen- 
tral source and in photoionization equilibrium. The value 
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Table 1. Journal of the BeppoSAX observation of NGC 3783. The different flux level states defined in Fig. |l| are 
indicated in each row. The spectral variability is discussed in Sect. ^. 



State 


Exposure (ks) 


"Count Rate (cts/s) 


f S/N 


•c 


*X 2 (dof) 


*P x2 


TOT 


153870 


0.758 ±0.002 


341 








LI 


27403 


0.586 ±0.005 


127 








HI 


35630 


0.823 ±0.005 


171 








Ml 


64257 


0.709 ±0.003 


213 








M2 


13311 


0.720 ±0.007 


99 








H2 


26791 


0.851 ±0.006 


151 








H = H1+H2 


62421 


0.835 ±0.004 


228 


1.02 


89.63(82) 


0.26 


L=M1+M2 


77568 


0.713 ±0.003 


235 


1.07 


109.5(82) 


0.02 


ratio H/L 








1.254 


222.3(82) 


<< 10" 3 



° In 1.5-10 keV. 

' Signal to noise ratio in 1.5-10 keV. 

* best fit value when the ratio of the two spectra is modeled with a constant 

* x 2 value for a constant model with const=C 

* Probability of exceeding x 2 ■ 
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Fig. 1. LECS, MECS and PDS lightcurves of BeppoSAX observation of NGC 3783. The selected different flux level 
states are also indicated. 



of the equivalent width for the OVII emission line pro- 
duced in a gas characterized by lg(U) — 1.29±q }| and 
lg(N H ) = 21.85±[J; 1 2 (as t hat observed by Bep poSAX), 



should not exceed ~ 10 eV ( Nicastro et al. 1999 ). If colli- 
sional ionization is also important, the emissivity from the 
gas strongly increases (Nicastro et al. 1999). In this case 



the blend of the OVII Ka triplet is more prominent and 
the predicted EW would increase up to 70 eV, still lower 
than measured. 

The HETG on board Chandra allowed to resolve the 



OVII triplet in NGC 3783 flKaspi et al. 20001) ; the total 
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BeppoSAX observation of NGC 3783. 
Model: Power Law + Cold Galactic Absorber 
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Fig. 2. LECS, MECS and PDS data (upper panel), and data/model ratio (lower panel) in the case of an emission 
continuum fitted with a power law absorbed by a cold galactic gas (model A). 



equivalent width was ~ 10 eV (even if poorly constrained) 
i.e. much smaller than that measured by BeppoSAX. 

All these considerations suggest that most of the low 
energy excess, which modeled with an emission line, was 
instead due to a soft X-ray continuum. We therefore fitted 
the excess at low energies with a black body component 
(model E). The best fit parameters are shown in Table 0. 
Comparing model E with model C (in which no oxygen 
line was included), we obtain that the black body compo- 
nent is required at the 99.99% confidence level (F=35.9). 

In the upper plot of the right panel of Figure ^| we 
show the ratio data/model E. When we add a gaussian 
emission line to the model E (model F) , with the intensity 
fixed at the Chandra value, we find that this component 
is required at only the 77% confidence level. A good fit 
(x 2 /dof = 105.1/102) is also obtained fitting the excess 
at low energies with a thermal bremsstrahlung (the ratio 
data/model for a thermal bremsstrahlung is shown in the 
last panel on the right in figure [j], model G), with tem- 
perature kT=0.46±0.06 keV. However, in this model the 
value of the relative reflection is smaller (i?=0.34 ±0.23) 
than expected from the iron line EWp e (~ 200 eV). We 
therefore conclude that the BeppoSAX spectrum is best 
fitted by model F. The best fit model F is shown in Figure 
H where each spectral component is labelled. We plot in 
Figure ^| the intervals of confidence for the parameter Ec 
and r (left panel), R and L (middle panel) and <jp e and 
Efc (right panel) 

4. Spectral variability analysis 



4-1- A model independent variability analysis 

To investigate spectral variability, we extracted spectra 
from different flux level states of the source during the 
BeppoSAX observation (see Table |l|, see Figure [l]) . 

The states we selected were: High 1 (HI), Low 1 (LI), 
Medium 1 (Ml), High 2 (H2) and Medium 2 (M2). To 
increase the statistics and to get comparable S/N ratios 
of the spectra extracted form the different flux levels, we 
added together those with similar flux. We fitted the ratios 
of the various spectra with a constant model, r(E) — C, 
to search for spectral variability. The best fit values for C 
are shown in Table [IJ The ratios H1/H2 and M1/M2 are 
both consistent with a constant value C~l (see Table Q). 
Then we create a new high state H=H1+H2, and a new 
low state L=M1+M2. In Figure ^ we present the ratio 
between H and L spectra. We can reject the hypothesis 
of constant value for this ratio at > 3er confidence level 
(see Table |l|). The plotted curves represent the expected 
spectral ratio in different cases. The variations can not be 
accounted for by a change of the warm absorber status 
(e.g. ionization parameter U). In fact we detected devia- 
tions of the spectral ratio from the constant model also 
above 2 keV (see Figure ^), where the warm gas is trans- 
parent to the radiation and no features can be produced 
by this material. In Figure ^ we also plotted the expected 
ratio if the only parameter which changes in the different 
flux states is the ionization of the warm gas (dash line). 
We suppose that the ionization structure is dominated by 
the intrinsic flux between 0.5 and 2 keV, this because the 
opacity of the gas is dominated by K shell absorption edges 
due to O VII and O VIII (the two most abundant ions of 
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Fig. 3. Ratio data model for the different spectral emission models employed to fit the NGC 3783 BeppoSAX data. 
See Table 0. 



the most abundant element). So we rescaled U assuming 

Uh = UL{Lg 5 _ 2keV I Lq 5 _2kev) ■ In tnis case tne absolute 
normalization of the Compton reflection is kept constant 
in the H and L states. The other curves in Figure || rep- 
resent the expected ratio if the observed variations were 
produced by: 

1) a change of the intrinsic slope of AT = 0.1 (dot 
- dash curve). The intensity of the reflected radiation is 
assumed to remain constant (i.e. the relative reflection R 
is variable). 

2) a change of the intrinsic slope as well as the re- 
flected component (solid line). In this case the reflection 
component is assumed to promptly respond to the inten- 
sity variation of the continuum (i.e. relative reflection R 
is constant). 

Looking at Figure |[ it is clear that the limited statis- 
tics does not allow us to distinguish between the two cases 
of constant or variable reflection. The observed behaviour 
can be reproduced by a variation of the intrinsic slope. 
This will be confirmed with a spectral analysis on the 
separate flux states, discussed in the next section. 

4-2. Spectral variability: the different flux level spectra 

We performed a spectral analysis of the different flux spec- 
tra, extracted as shown in Figure [l] and in Table |l|. The 



best fit model employed in this analysis is the one giving 
the best fit to the integrated spectra, i.e. the model F. 

We confirm that the observed variations can be ex- 
plained by a change of the intrinsic continuum. In partic- 
ular the source steepens (Ar ~ 0.1, see Figure Table 
||) when brightens. The limited statistics in the separate 
flux states did not allow us to detect variations of E c (see 
contour plot in Figure 0). No variations in the ionization 
status of the warm gas was observed. 



5. Discussion 

5.1. Intrinsic Continuum and spectral variability 

The presence of a high energy cut-off was required to 
fit the X-ray/7-ray spectra of a sample of Seyfert ob- 
served with different satellites (Zdziarski et al. 1996). The 



E c = 340±J|$ keV detected in NGC 3783 by BeppoSAX 
is consistent with the values observed in other Sy Is of 
the BeppoSAX sample (ranging between 70 keV up to 
300 keV, Perola et al. in preparation). 

The spectral variability detected in our observation 
can be accounted for by variations of the intrinsic slope. In 
particular the intrinsic spectrum steepens when the source 
brightens (AT - 0.1 with (AF/F) 2 -iokeV ~ 20%). This 
behaviour was observed also in other sources: NGC 5548 
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Table 2. NGC 3783: Spectral emission models 







Moel B 


Model C 


Model D 


Model E 


Model F 


Warm 


lg(U) 




i.eo±8:8J 


i.29±8:« 


i.29±r 9 


1.33 ±0.10 


Absorber 




OO oo i 0-03 


22.31±8:8e 


21.85±8:i4 


22.28±8:ol 


22.26±8:8I 


Intrinsic 


E c (keV) 




299±fgl 


143±H 


307±M° 


340±?o? 


Continuum 


r 


1.90±0.04 


i.89±8:S 


i.69±8:8S 


1.86±0.03 


i-87±8:8? 




°Norm. 


2.2±0.1 


2.1±0.1 


1.6±0.1 


2.1±0.1 


2.1±0.1 


Reflection 


*R 


o.75±8:?? 


0.95±g;fl 


0.31±8;?^ 


0.71±g;i 


o.ro=b8:to 


and 


E L (keV) 


6.41 ±0.04 


6.41 ±0.04 


6.40 ±0.04 


6.39 ±0.09 


6.39 ± 0.09 


Iron Line 


gl (keV) 


o.53±8:?i 


o.49±EJ:i 


o.39±8:? 


0.38±oiog 


o.36±8:81 




EWl (eV) 


300 ±45 


270 ±60 


216 ±50 


210±45 


190±47 


OVII 


E (keV) 






0.51 ±0.02 


_ 


0.52±0.06 


Line 


(JO (keV) 


- 


- 


o.i3±8:8§ 


- 


< 0.2 




EWo (eV) 






240 ±27 




10 (frozen) 


Soft 


kT (keV) 








0.21 ±0.01 


0.22±0.03 


Emission 


t1 "Norm. 








4.3± ) :1 


3.4±oJ 




X 2 /dof 


189.6/108 


176.2/107 


111.2/104 


104.7/105 


101.7/103 




'F value 


59.7 (99.9%) 


8.14 (99.5%) 


20.3(99.9%) 


35.9(99.99%) 


1.5(77%) 



* In 10~ 2 photons keV -1 cm~ 2 s' 1 , a 1 keV. 

* R = Af2/27r is the Relative reflection, the solid angle subtended by the reflecting medium. The inclination angle between 
the line of sight and the reflector is assumed i=30°. 

T' In 10 -4 Z/3g/D 2 , where L39 is the source luminosity in units of 10 39 ergs/sec and D10 is the distance to the source in units 
of 10 kpc. 

t In respect to model A (see text) for model B, in respect to model B for model C, in respect to model C for model D, in 
respect to model C for model E, in respect to model E for the model F. 



( [Nicastro et al. 2000[ ), NGC 4151 (Piro 2002), NGC 7469 
(Nandra et al. 2000), MCG-6-30-15 (Vaughan & Edelson 
2001), and IC 4329A (Done et al. 2000). 

A two-phase model involving a hot corona emit- 
ting medium-hard X-rays by Comptonization and a cold 
optically thick layer which provides the soft photons 
to be Comptonized (Haardt & Maraschi 1991, Haardt 
et al. 1997), can explain the observed intrinsic variability. 
In this model the physical parameters of the corona can 
drive the spectral variations observed in the hard-medium 
energy range, while those in the soft band are driven by 
the temperature Tbb of the soft photons emitted by the 
cold layer. 



If the corona is not pair dominated, i.e. the compact- 
ness parameter l c — c £^ c a^ < 10, the two phases model 
predicts a noticeable spectral variability even when the 
luminosity variation is less than a factor two. Our vari- 
ability study suggests that this is the case for NGC 3783. 
We measure Ar ~ 0.1 (see Figure (?] and Table |^), with no 
variation of the Comptonization luminosity (0.1-200 keV) 
in the H and L states: ~ ~ 10 44 erg/s. We can 
put a lower limit to the compactness parameter. The time 
scale over which the flux changes by a factor 2 is greater 
than 5 days, so we obtain l c > 0.2. A detailed Comp- 
tonization code was employed to fit the data of the Bep- 
poSAX observation of NGC 5548 ( [Petrucci et al. 20001) . 
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The BeppoSAX broad-band spectrum and variability of the Seyfert 1 NGC 3783 
NGC 3783: Best fit spectrum: Model F. 




channel energy (keV) 

Fig. 4. Unfolded spectrum and the best fit model F (See Table ||. The contributions to the model of the various 
additive components are also plotted 

NGC 3793: Total spectrum NGC 3783: Total spectrum NGC 3783. Iron Line 




Photon Index Photon Index ^ M 

Fig. 5. la, 2(7 and 3(7 confidence levels E c vs T (left panel), R vs T (middle panel) and crpe vs Ei? e (right panel) in 
the total spectrum. 



The spectral variability observed in that source could by 
accounted for by a change of the spectral slope AF ~ 0.2 
with small hard luminosity changes. The spectral soften- 
ing in the high state is explained (when a realistic model 
of Comptonization is employed) by a decrease of kT e due 
to an increase of the soft photon flux. The anticorrelation 
between the slope and the high energy cut-off is not re- 



vealed when a cut-off power law is used as a zero order 
approximation of a Comptonized spectrum. Also in the 
case of NGC 5548 the trend of variability suggests a hot 
plasma which is not pair-dominated (and the estimated 
lower limit for the compactness was l c > 0.2). In the case 
of NGC 3783 an accurate Comptonization model was em- 
ployed only to fit the total spectrum, with no spectral 
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NGC 3783: Spectral ratio H/I, 




1 10 

Energy (in keV) 



Fig. 6. Ratio of spectra extracted by different flux level states. The plotted curves represent the expected behaviour 
of the ratio when the variations are produced by change of the intrinsic continuum, either with relative reflection 
R=const (solid line) or with relative reflection R variable (dot-dash curve), and by a change of the warm gas (dash 
curve). The dot line is the best fit value when the ratio is fitted with a constant. See text for detail. 
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Fig. 7. NGC 3783. la and 2a confidence levels E c vs T in the high (H) and low (L) states. 



variability analysis ( |Petrucci et al. 200f[ ). The estimated 
coronal parameters are substantially different between the 
two models (the coronal temperature kT e = 265±5 5 keV 
is higher than that measured with our simpe cut-off power 
law model E c /2). This is due to the presence of an energy 
break in the Comptonization model, produced when the 
anisotropy of the soft photon field is taken into account. 



5.2. Warm absorber and soft emission 

The low energy spectrum of NGC 3783 is character- 
ized by the presence of the warm absorber. The com- 
plex structure of this gas was inv estigated in the Chan- 
dra observation ( Kaspi et al. 200f ). The model consisted 
of two absorption/emission components with the same 
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Table 3. NGC 3783: Spectral variability, best fit parameters in the different flux states in the model F. 







Low State: L 


High State: H 


Warm 


lg(U) 


i.28±8:o? 


i.38±8:w 


Absorber 




22.27±0.11 


22.32±8:!o 


Intrinsic 


E c (keV) 


> 250 


300±?g!] 


Continuum 


r 


1.83±0.04 


i.93±8:8S 




°Norm. 


1.9 


2.5 


Reflection 


°R 


0.58±8;il 


o.77±8:£ 


and 


E L (keV) 


6.39±0.07 


6.39±0.06 


Iron Line 


gl (keV) 


0.37±S];? ( 7 


0.25±8;i 




EWl (eV) 


204±78 


163±70 


OVII 


E 2 (keV) 


0.5±0.2 


0.5±0.2 


Line 


o-2 (keV) 


(frozen) 


(frozen) 




EWo (eV) 


34±51 


30±41 


Soft 


kT (keV) 


o.2i±8;8i 


0.22±8;il 


Emission 


*Norm. 


4.5±J:? 


4.6± :| 




x 2 /doi 


85.26/102 


73.20/102 



* In 10~ 2 photons keV" 1 cm" 2 s _ \ a 1 keV. 
° with cos(i)=30° 

* In 10 _4 Z/39/-Dio, where L39 is the source luminosity in units of 10 39 ergs/sec and D10 is the distance to the source in units of 
10 kpc. 



column density lg(iVff) = 22.1, but different ionization 
parameter: lg(U oxvgen ) = —0.745, with a covering fac- 
tor f=0.3 and \g{U oxygen ) = —1.745, with a covering 
factor f=0.5 (these values are computed between 0.538 
keV and 10 keV). We fitted the BeppoSAX data with 
a grid of single-zone warm absorber models build with 
CLOUDY QFcrland 1997| ). These models include only the 
transmitted spectra, and do not include the contributions 
of both gas emission and resonant absorption. The ion- 
ization parameter, computed from E > 1 Ryd, in the 
total spectrum was lg(J7) = 1.33 (see Table g). Taking 
into account the ionizing continuum used to fit Chandra 
data, we get the relation U oxyflen =0.114 U(0.1-10 keV) 



(Kaspi et al. 2001); in the case of BeppoSAX spectrum 
we employed the ionizing continuum from Mathews & Fer- 
land (1995) (listed in Table J|) to rescale the value of U(> 
1 Ryd) in the 0.1-10 keV energy range. In this case we 
found U(> 1 Ryd) - 50 U(0.1-10 keV) (see also George 



et al. 1998a). The ionization parameter during our obser- 
vation is marginally consistent with that of the low ioniza- 
tion gas observed in the Chandra spectrum t/sAx(0.1-10 
keV)=0.4±0.1. No variation of the low-ionization compo- 
nent was observed between BeppoSAX and Chandra ob- 
servations, not surprisingly as no flux changes in 0.5-2 keV 
(F .5_2fc e y — 2 x 10~ n erg cm' 2 s~ x ) between the two 
observations have been detected. 

The high ionized absorption features of the heavier el- 
ements produced in the more ionized component are too 
small to be observed given the limited LECS sensitivity 
and energy resolution. The warm absorber parameters ob- 
tained with BeppoSAX are derived mainly from the deep 
OVII and OVIII absorption edges produced in the low 
ionization gas. The OVII triplet observed by Chandra at 
~ 0.5 — 0.6 keV is only marginally required by our data 
(Pi?=77%). This feature was also observed in the Bep- 
poSAX spectrum of NGC 5548 ( |Nicastro et al. 2000|) , and 
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later on confirmed by Chandra (Kaastra et al. 2000). An 
additional absorption line (at ~ 1 keV) corresponding to 
the blend of many iron L-shell absorption lines observed 
by Chandra, is very marginally required by the model 
(Pir=52%). When we add this component the warm ab- 
sorber parameters as well as those of the intrinsic contin- 
uum do not change. 

A reanalysis of the ASCA observation of 1996 taking 
into account (1) the inclusion of the reflection compo- 
nent as measured by the RXTE observation simultane- 
ous with those Chandra in 2000 and (2) the use of the 
Yaqoob (2000) method to account for the degradation of 
the SIS detectors, was performed by Kaspi et al. (2001). 
The intrinsic continuum they found is fully consistent with 
that measured by BeppoSAX (r=1.85±0.03, A(l keV) 
=2.21xl0~ 2 ph cm~ 2 s" 1 keV -1 ), and the warm absorber 
is consistent with the low ionization component observed 
by Chandra (J7 X =0.133±0.003, lg Afa=22.15±0.03). 

A significant soft emission well modeled by a black 
body component with a temperature kT=220±30 eV was 
detected in the BeppoSAX observation. No variation of 
the soft component emission was detected. We may as- 
sociate this emission to the accretion disc around the 
central source. The IUE monitoring of NGC 3783 from 
1978 to 1985 detected a UV-bump from 5150 A up to 
the soft X-rays. The observed bump could be attributed 
to a thermal emission from an optically thick medium. 
When the UV bump was reproduced with a sum of black- 
body spectra, the highest disc temperature was 280,000 K 
( [Alloin et al. 1995| ). 



Table 4. Ionized continuum from Mathews & Ferland 
(1987) employed to build a grid of single-zone warm ab- 
sorber model to fit the low energy spectrum of NGC 3783 



v(Ryd) 




slope a 


0.206 


2.6576 


0.5 


1.743 


2.194 


1 


4.130 


1.819 


3 


26.84 


-0.6192 


0.7 



5.3. Compton reflection and Iron line 

The broad band capability of BeppoSAX allowed us to 
detect the Compton reflection together with a broad Fe 
Ka emission line (Figure ||). The energy of the iron line 
E Fe = 6.39 ± 0.09 keV, is consistent with Fel-XII. In the 
model F (see Table this features is well reproduced 
by a broad gaussian component with op e = 0.36±Q;[]g 
keV and EWfc=190±47 eV. A broad and asymmetric line 
is expected if the X-ray reflection is occurring in the in- 
ner region of an accretion disc (Fabian et al. 2000, Matt 
et al. 1992). However there is evidence that in NGC 3783 
as well as in other Seyfert 1 galaxies, the line is actually 



composed of two components, one broad and variable and 
the other narrow and constant. To support this idea we 
recall that in the Chandra observation on NGC 3783 a nar- 
row unresolved iron line (with EW=115±36 eV) was de- 
tected (Kaspi et al. 2001). The upper limit for the FWHM 
was 3250 km s , consistent with this line originating out- 
side the BLR. The line flux is consistent with predictions 
from models of torus emission. The Chandra data could 
not constrain any broad Fe Ka emission though they were 
consistent with the broa d line found by AS CA in this 
source (EW=230±70 eV, |George et al. 1998b| ). 

If we fit the iron line with a narrow component with 
the parameters fixed to those detected by Chandra (E„; = 
6.39 keV, I„; = 6.6 x 10 -5 ph cm~ 2 s^" 1 ), we find strong 
residuals around 6 and 7 keV (see Figure [|, left panel) . If 
we add a second broad gaussian component the residuals 
disappear (see right panel in Figure ^). The broad com- 
ponent is required at 99.9% c.l. (following a F-test adding 
three interesting parameters), with au = 0.72±Q-|y keV 
EW W = 115 ±76 and E hi = 6.27±g;§;5 keV, i.e. correspond- 
ing to neutral iron. These values are fully consistent with 
those of the ASCA observation in 1996 when the ASCA 
data were modeled assuming the underlying continuum 
derived by the analysis of a RXTE observation (simulta- 
neous with those of Chandra in 2000, Kaspi et al. 2001), 
to constrain the reflection component. 

The simultaneous BeppoSAX/XMM observation of 
the NLSyl NGC 5506 QMatt et al. 2001b[) showed that in 



that case the bulk of the Compton reflection has to be 
associated with the narrow iron line, as in the the case of 
NGC 4051 (Guianazzi et al. 1998). In other objects ob- 
served by XMM (Mkn 509, Pounds et al. 2001; Mkn 205, 
Reeves et al. 2001), evidence for ionized accretion disc 
was also found. An ionized disc reflection model was also 
employed to fit the data of NGC 3783. This source was 
included in a small sample of Sy Is observed by Bep- 
poSAX analysed with the ionized disc refle ction model of 
Ross & Fabian (1993) ( |De Rosa et al. 2002j ). In this model 
the most important quantity in determining the shape 
of the reflected continuum is the ionization parameter 
£ = 4:TrF x /nff, where F^, is the X-ray flux (between 0.01- 
100 keV) illuminating a slab of gas with solar abundances 
and constant hydrogen number density n#=10 15 cm~ 3 . 
The Ka iron emission line as well as the emission features 
from the disc at low energy, are included in the model. In 
the fit, we applied to the spectrum a relativistic blurring 
appropriate for a Schwarszchild geometry assuming a disc 
emissivity law (Fabian et al. 1989). During the fit a nar- 
row iron line component with all the parameters linked to 
those of Chandra observation was included. This analy- 
sis showed that the reflected spectrum of NGC 3783 can 
be produced (x 2 /dof=104/101) by a highly ionized disc 
characterized by: inner radius (pegged to the lower value) 



=6 r , outer radius r Q 



7+ 1 



Yg With Tg- 



=GM/c 2 the 



gravitational radius, ionization parameter lg£ = 2.73± 
and inclination angle i = 35.6±22 degrees (see Figure |l0|) 



0.04 

'5 
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NGC 3783. Iron Line & Compton Reflection 




10 20 
channel energy (keV) 

Fig. 8. Ratio data/best fit model when the normalization of the iron line and the Compton reflection component are 
set equal zero. 



NGC 3783. Narrow iron line. Residuals 



channel energy (keV) 



NGC 3783. Narrow & Broad iron line. Residuals 



channel energy (keV) 



Fig. 9. Left panel. Residuals to narrow iron line when all the line parameters fixed to those of Chandra observation. 
Right panel. A similar plot to left panel with the addition of a broad gaussian component. 



In this case the observed soft X-rays excess can be ac- 
counted for by the disc emissivity and no additional com- 
ponent at E < 2 keV is required. This suggests that the 
Compton reflection observed in NGC 3783 can be com- 
pletely associated with the broad Fe component suppos- 
ing that both were created in a relativistic disc. Unfortu- 
nately, the energy resolution and sensitivity of BeppoSAX 
are not adequate to disentangle both the Fe line compo- 
nents and simultaneously the fraction of reflection associ- 
ated with them. A simultaneous BeppoSAX/XMM obser- 
vation would be the best way to achieve this purpose. 



6. Summary 



We have analysed the long BeppoSAX observation of the 
Seyfert 1 galaxy NGC 3783. We detect a high energy cut- 
off of the intrinsic power; the value of E c —340±\q^ keV 
is consistent with those observed so far in a sample of 
Seyfert galaxies, which ranges from 70 keV up to 300 keV. 
Spectral variability was clearly present, well modeled by 
a steepening of the intrinsic continuum (AT ~ 0.1) with 
no total luminosity change. This behaviour is consistent 
with a two-phase disc-corona model in which with the hot 
plasma is not pairs dominated (as observed also in the 
Seyfert 1 galaxy NGC 5548). 
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Fig. 10. The unfolded BeppoSAX/MECS spectra (EF E ) of NGC 3783 is plotted when the ionized disc reflection model 
of Ross & Fabian (1993) is employed to fit the data. A narrow iron line component with all the parameters linked to 
those observed by Chandra is added to the model. A broad disc component is also clearly present. 



The presence of a warm gas is confirmed by our 
analysis. Its ionization status is consistent with the low- 
ionization component detected by Chandra. NGC 3783 is 
one of the few sources in which BeppoSAX detected a soft 
excess. We reproduced the soft emission with a black body 
component with a temperature kT^0.2 keV. 

The Compton reflection hump and the iron emission 
line are clearly detected. In addition to the narrow un- 
resolved iron component observed by Chandra, a broad 
iron line component (a = 0.72±J'27 keV) is required to fit 
the BeppoSAX data. The observed reprocessing features 
can be reproduced with a highly ionized accretion disc 
(lg£ = 2.73±q;[]|). In this case the soft excess observed at 
E < 2 keV can be accounted for by the disc emissivity and 
no additional soft X-rays component is required to fit the 
spectrum. 
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